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Abstract. A new form of Bethe ansatz equations is introduced. A version of a sepa- 
ration of variables for the quantum 5(3 Gaudin model is presented. 



1. Introduction 

The separation of variables for the quantum SI2 Gaudin model was constructed by 
Sklyanin in [Sk~j . In this paper, we give an analogue of Sklyanin's construction for the 
Lie algebra sl 3 . 

We were inspired by Stoyanovsky's paper |5tj . in which the author uses Sklyanin's 
change of variables to establish a relation between the Knizhnik-Zamolodchikov equa- 
tions |KZj and the Belavin-Polyakov-Zamolodchikov equations |BPZj . and to construct 
integral formulae for solutions to the BPZ equations. 

1.1. The paper is organized as follows. In Section |21 we introduce notations and recall 
the definition of the KZ equations and Gaudin model. 
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1.2. In Section El we recall the definition of the master function and canonical weight 
function. In Theorem 13.4.11 we recall the main fact of the Bethe ansatz method: if t is a 
critical point of the master function, then the value at t of the canonical weight function 
is an eigenvector of the Gaudin Hamiltonians, see |RVj . cf. |Balt IB"a2] . 

The critical point equations for the master function are called also the Bethe ansatz 
equations for the Gaudin model. 

1.3. In Section we discuss different forms of the Bethe ansatz equations for the Lie 
algebras 5^ and SI3. 

For sl 2 , there are two ways to describe solutions to the Bethe ansatz equations. The 
original way: a solution is a collection of numbers t = (t^\ . . . , t[ ), satisfying the critical 
point equations (JT2J). The second way: given the numbers zi, . . . , z n , (Aj, oti), (Aj, Aj), a 
solution is a polynomial P(x) = n!Li(4 — x) of one variable and numbers /ii, . . . 
/ii + • • ■ + yU„ = 0, satisfying the differential equation 



(Aj, A,-) . 

v - 3) =0 



i=i i=i jjti ■> 

The equivalence of two descriptions is a classical fact which goes back to Stieltjes, see 
StiJ and Sec. 6.8 in [5z] . 

For sis, there are four ways to describe solutions to the Bethe ansatz equations. In 
this introduction, we mention only two of the four: the original way and the new way 
suggested in this paper. 

The original way: a solution is a collection of numbers t = (t^\ . . . , tf^\ tf\ . . . , tf^), 
satisfying the critical point equations (|15|) . The new way: given the numbers z\, . . . , z n , 
(Aj, «i), (Aj, a 2 ), (Aj, Aj), a solution is a pair 

Pi{x) = U(tf ] -x) , P 2 (x) = flitf-x) 
1=1 1=1 

of polynomials of one variable and a set of numbers /ii, . . . ,/x n , (J>i + • • • + (J> n = 0, 
satisfying the differential equation 

(2) p['p 2 - p[p< 2 + Pl p» - p[p 2 £ { -^- - p i p 2 E + 



1=1 1 1=1 1 



n 



t f /• y . t J y . y . 

■ -1 ^ ■ , ■ ~% ^7 

1=1 J^l J 

see Theorem 14.2.31 This theorem is the first main result of our paper. 



1.4. In Section we describe Sklyanin's separation of variables for the sl 2 Gaudin 
model, following the exposition in |5t] . 
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The general problem in the Gaudin model is to find eigenfunctions of the commuting 
Gaudin Hamiltonians Hi(z), i = l,...,n, where Hi(z) + ■ • ■ + H n (z) = 0. In suit- 
able coordinates, the Hamiltonians are differential operators acting on polynomials in n 
variables , . . . , x^ n ' . 

The eigenfunction equations are 

Hi{z)F^\...,x^) = ^F(x^,...,x^) , z = l,...,n, 

where the eigenvalues satisfy the equation fix + • • • + \i n — 0. The famous Sklyanin's 
change of variables (fT%J) from variables x^\ . . . , x^ n > to new variables u, y^\ . . . , y^ 1 ^ 
transforms the eigenfunction equations to the following equations for the unknown poly- 
nomial F(x(u, y)), 



i=l y 4 i=l y 1 fe^j 4 fc 

j = 1, . . . , n — 1, see |Sk| 1ST] and Theorem 15.3.11 

This Sklyanin's statement has three interesting features. 

The first is that the variables have separated: the j-th equation depends on variable 
only and does not depend on u and other variables y^\ 

The second feature is that the differential operator is the same in all equations. 
Namely, the differential operator for an index j' can be obtained from the differential 
operator for the index j by replacing y^> with y^>. 

The third feature is that the differential operator in Sklyanin's equations is the same 
as the differential operator in the Bethe ansatz equation 

Sklyanin and Stoyanovsky also consider the canonical weight function. In Sklyanin's 
variables, the canonical weight function takes the form 

V(t,z,u,y) = u h ^-^ ^ ' 



where l\ is a nonnegative integer and P(x) = Yli=i{b ' ~ ' 



P(z l )...P{z T) 

The canonical weight function, as a function of u, y^\ . . . , y^ 1 ', is the product of 
functions of one variable. This is another manifestation of separation of variables. 

Now if one looks for a value of the parameters t such that z, u, y), as a function of 
u, y, becomes an eigenfunction of the Gaudin operators, then one gets the Bethe ansatz 
equation for the unknown polynomial P(x). 

Hence, if for some t, the function \l/(t, z, u, y) is an eigenfunction, then t satisfies the 
Bethe ansatz equations. 

1.5. In Section IH1 we describe an analog of Sklyanin's statements for sl 3 . 

Again the general problem is to find eigenfunctions of the Gaudin Hamiltonians 
Hi(z), i — 1, . . . , n, where Hi(z) + ■ • • + H n (z) = 0. Now the Hamiltonians are differ- 
ential operators acting on polynomials in 3n variables Xi,X2~\x^\ . . . , x^ ,x^\x^ . 
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The eigenfunction equations are 

Hi(z)F(x) = HiF(x) , i = l,...,n, 

where the eigenvalues satisfy the equation /ii + ■ ■ • + \i n = 0. 

We make a change from variables x^\x^\x^\ , n , £3 to new variables 

u\, w 2 , M3, Ui'\y { 2\y { i\ ■ ■ ■ , Hi 1 , V2 , 2/3™ , see formula (|27ji analogous to Sklyanin's 
formula. 

We define the degree of a polynomial in C[u, y] as its degree with respect to variable 
u 3 . 

On the affine space with coordinates u, y, we consider the affine subspace 

D = {(u,y) I = j = 1, . . . , n - 1} . 

It turns out that the Gaudin Hamiltonians have the decomposition 

H l (z)= HiMa + Hiiz^o + HiizUo, 

where the operator Hi(z) preserves the degree, the operator Hi(z) >0 increases the 
degree by one, and the operator H^z)^ decreases the degree by one. Moreover, it turns 
out that for any F e C[u, y], the restriction to D of the polynomial Hi(z) <0 F is zero. 

These remarks show that the eigenfunction equations become "upper triangular" with 
respect to the degree decomposition after restriction to D. 

Let F(u, y) be an unknown eigenfunction with eigenvalues /ii, . . . , \i n . Let F(u, y) = 
F(u,y)o + F{u,y)i + . . . be the degree decomposition. According to our change of 
variables, 

where l\, I2 are nonnegative integers and / is a polynomial depending on variables 

Vi\v2\ 3 = !, • • -,n- 1, only. 

The triangularity remark implies the following equations for the degree zero part of 
the eigenfunction, 

(4) ((Hi(z) -fH)F ) =0, z = l,...,n. 

1.5.1. Theorem. After suitable renormalization, equations take the form 

(K\ ( ° 2f ° 2f ° 2f 9f V (A " Ql) 

^ > \ a..U)a..U) (?) 



WW ^ J 2/i JJ - * 

g/ (Aj,a 2 ) , 1 / \- (A^A fc ) 

~ 0) 2^ u) _ + J 2^ u) _ z . [ + 2^ z _ z 

°y2 i=i U\ z i i=i tfi z * fc^i * K 



= . 



This is our second main result, see Theorem 16.5.41 

As Sklyanin's equations © for sl 2 , our equations (J3j) have three interesting properties. 
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The first is that the variables have separated: the j-th equation depends on variables 
Vi\ onr y an d does not depend on u and other variables y± \ y% ^ 

The second property is that the differential operator is the same in all equations. 
Namely, the differential operator for an index j' can be obtained from the differential 
operator for the index j by replacing 2/1,3/2 with y± \ y!f \ respectively. 

The third property is that the differential operator in equations © is the same as the 
differential operator in the Bethe ansatz equation (J2J). 

Then we consider the canonical weight function z, u, y) and its weight decompo- 
sition = \J/ + ^1 + It turns our that the degree zero term has the form 



(6) V(t,z,u,y) Q = u^u 



X) 



where I1J2 are nonnegative integers, P\{x) = Yli=i(H x )> P2\ x ) = lLLil^i 
Moreover, it turns out that the higher degree terms of \l/ are uniquely determined by 
the degree zero term \l/ , see Theorem 16.6.31 and Section IB. 7. 31 

According to formula (jUJ), the degree zero term \l/o of the canonical weight function is 
the product of functions of one variable. This is another manifestation of separation of 
variables. 

Now assume that one looks for a value of parameters t such that ty(t, z,u,y), as a 
function of u, y, becomes an eigenfunction of the Gaudin operators. Then, by Theorem 
11.5.11 one gets the Bethe ansatz equation (j2J) for the unknown polynomials -Pi(x), P2{x). 

Hence, if for some t, the function z, u, y) is an eigenfunction, then t satisfies the 
Bethe ansatz equations. The precise statement see in Theorem 16.7.41 



1.6. In the appendix we consider the Bethe ansatz equations for the Gaudin model 
associated with an arbitrary Kac-Moody algebra of rank r. We introduce a polylinear 
differential equation for a collection Pi(x), . . . , P r {x) of polynomials of one variable and 
a collection of numbers . . . , fi n , jii + • — h /i„ =0. That equation is an analog of the 
differential equation (j2J). Under certain conditions, we show that the union of roots of 
polynomials P±(x), . . . , P r (x) form a solution to the Bethe ansatz equations if and only 
if Pi (x), ... , P r (x) and /xi, . . . , /i n satisfy that polylinear differential equation. 

This fact may be considered as a generalization of Stieltjes' Lemma 14.1.21 (see also 
|Stij . and Sec. 6.8 in jSij) to an arbitrary Kac-Moody algebra. 



1.7. In the next paper we plan to extend the results of this paper from 5(3 to other Lie 
algebras. 



1.8. This work has been started when the fourth author visited Universite Paul Sabatier 
in Toulouse in May- June of 2006. He thanks the university for hospitality. 
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2. KZ EQUATIONS AND THE GAUDIN MODEL 

2.1. Lie algebra sl r+ i. Consider the Lie algebra gl r+1 with standard generators e at b, 
a,b — 1, . . . , r + 1. Set h a = e a ^ a — e a+ i ja+ i for a = 1, . . . , r. Consider the Lie subalgebra 

— Slr+i 

generated by elements e 0j &, a ^ b, and h a , a — 1, . . . ,r. Let g = n © f) © n + , 

n_ = © a >6 C • e a , b , fj = ©^ =1 C • /i a , n+ = © a<6 C • e a , b , 

be the Cartan decomposition. Let a a e f) *, a = 1, . . . , r, be simple roots. 

Fix the invariant scalar product on g such that (/ij, hi) = 2. The scalar product 
identifies g and its dual g*. 

2.1.1. The Casimir element f2 e g® 2 is the element © Xj, where (x^ is an or- 
thonormal basis of g . 

2.1.2. For = sl 2 , 

1 

Q = e 2 ,i © ei >2 + ei )2 © e 2 ,i + -h © /ii . 

For g = sl 3 , 

^ = e 2 ,i © ei )2 + ei )2 © e 2 ,i + 

e3,2 © e 2 , 3 + e 2 , 3 © e 3 ,2 + e 3 ,i © ei j3 + e lj3 © e 3) i + 

/ 2 1 \ / 1 2 \ 

/ii ® g/ii + ^h 2 J + h 2 ®{-h 1 + -h 2 J . 

2.1.3. Let 2/ (ti_) be the universal enveloping algebra of n_, 

C/(n_) = © IeZ r o C/(n_)[/] , 
where for I — . . . , / r ), the space C/(n_)[Z] consists of elements / such that 

r 

[f,h\ = (h,j2^)f- 
i=i 

The element f| i e ai) 6 j with a* > 6$ belongs to the graded subspace U(n~) [I], where 
I = £\ i W with 

J (i) = (0,0,...,0,l b ,l b+1 ,...,l a _!,0,0,...,0) . 

Fix an order on the set of elements e a ^ with r + l>a>b>l. Set e a ^ < e a > ^ if 
b < b' or 6 = b' and a < a'. 
For example e 2 ,i < e 3; i < e 3i2 . 

Then the ordered products ria>6 e afe b f° rm a graded PBW basis of t/(n_). 
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2.1.4. The grading of U(n—) induces the grading of C/(n_)® n for any n, 

f/(n_) 0n = ©, l/(n_)® n [l] , 

where U(n.)^ n [l] = ® l w + ... +l M =l U(n-)[l {1) }®- ■ -<g)C/(n_)[Z (n) ]. The PBW basis of tensor 
factors induces a graded PBW basis of U(xi-)® n . 

2.1.5. For A G f)*, denote by Ma the sl r+ i Verma module with highest weight A. 
Denote by v\ G M A its highest weight vector. 

For A = (Ai, . . . , A n ), A s G f) *, denote 

M A = M Al (8) • • • (8) M An . 
We have the weight decomposition 

M A = ©i eZ . o M A [Z] , 

where M A [(/i, . . . , / r )] is the subspace of vectors of weight Yl™=i ^-s — YH=i ^ a i- 
The PBW basis of U(n^)® n induces a graded PBW basis of M A . 
Let 

Sing M A [l] C M A [l] 

denote the subspace of singular vectors, i.e. the subspace of vectors annihilated by n + . 

2.2. KZ equations. The KZ equations on an M A -valued function I(z±, . . . , z n ) of com- 
plex variables z = (zi, . . . , z n ) is the following system of differential equations 

dl x ^ fiW) r 
k j— = — — — J, i = l,. ..,n, 

where k G C* is a parameter of the equations. 

The KZ equations will not be used in this paper, but it is useful to keep them in mind 
when discussing the Gaudin model. 

2.3. Gaudin model. Fix pairwise distinct complex numbers z±, . . . , z n . Denote 

H i( z ) = 2^ — ~ — ' i = l,...,n. 

These are linear operators on M A called the Gaudin Hamiltonians. 

The Gaudin Hamiltonians commute, [Hi(z), Hj(z)} = for % ^ j. 

The Gaudin Hamiltonians commute with the action of sl r +i on M A and hence they 
preserve the subspaces Sing M A [Z]. 

We also have Hi(z) H h H n (z) = 0. 
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3. Master function, canonical weight function, and the Bethe ansatz 

3.1. Master function, |S Vj . Let A = (Ai, . . . , A n ), A s E f) *, be a collection of sl r+ i- 
weights and / = (Zi, . . . , l r ) a collection of nonnegative integers. Set I = l\ + • • • + l r . 
Introduce a function of n variables z = (z\, . . . , z n ) and I variables 

t-(t {1) t {1) t (r) t {r) ) 
L — \ L i ■> ■ ■ ■ j ^ j ■ ■ ■ j L i i - ■ ■ i L i r ) 

by the formula 

(7) *(*;*; A; I) = \{ ( Zi - Zj )^ x 

l<i<j'<n 

r li n r 1 — 1 li 'i+i 

nniKf -^) ,A " n n (f-wnniM'-'rv. 

i=l j=l s=l i=l l<j<s<k i=l j=l k=l 

The function $ is a (multi- valued) function of t, depending on parameters z, A. The 
function $ is called the master function. 

3.1.1. The product of symmetric groups 

Si = S h x---x S lr 

acts on variables t by permuting the variables with the same upper index. The master 
function is ^-invariant. 

3.1.2. A point t with complex coordinates will be called a critical point of $( • ; z; A; I) 
if the following system of I equations is satisfied 

4> 1 — j (t;z;A;l) = , i = 1, . . . , r, j = 1, . . J; . 

3.1.3. Equations (jHJ) can be reformulated as the system of equations 



(A s , ct\) sr~^ 2 x 1 

( ) 1^ _ : ~ 1^ .a) _ M) + 2s f (i) _ i(2) = ' 

s=l L j z s s= i j s= i 

Ai) 2^ m _ Ai) + 2^ .(o _ .(i-i) + 2> _ u ' 

s=l L j z s s= i j SJ tj t-s s=1 t-j l s s=1 Lj L s 

71 / » \ / r lr — 1 ^ 

(A B ,a r ) 2 1 n 

2^ ,(r) 2^ ,(r) _ (r) + 2^ ,(r) _ (r-1) U ' 

s=l s s=l, s^j ^ ts s=l s 

where j = 1, . . . , l\ in the first group of equations, i = 2, . . . , r — 1 and j = 1, . . . , Zj in 
the second group of equations, j = 1, . . . , l r in the last group of equations. 
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3.2. Canonical weight function, jSVl IKSVj . The PBW-basis of C/(n_), defined in 
Section l2.1.3l is 

1 r io ~~ e 2,l e r+l,r J 5 

where Jo = {i a ,b}a>b runs over all sequences of nonnegative integers. 

Let Ffa . . . , Fj n be elements of the standard PBW-basis, Ij = {i 3 2 i, ■ ■ ■ , i 3 r +i r }- Set 
I — (Ii, . . . , I n ). The corresponding basis vector 

F/u = F/^m ® • • • ® F /n+1 t; An 

lies in M A [Z] if 

n i r+1 

( 10 ) EE E <6 = fora11 * = !,■•■,»• . 

j=l 6=1 a=i+l 

Denote by P(l,n) the set of all indices I corresponding to basis vectors in Ma [2]. The 
set {Fjv}i e pu n \ forms a basis of M A [Z]. 

3.2.1. For I G P(l,n), define the set 

S(I) = { (j, a,b,q) | l<j< n, l<6<a<r + l, 1 < 9 < < 6 } • 

For z = 1, . . . , r, define the subset 

S l (I) = { s = (j,a,b,q) G S(I) \ b<t<a} . 

Condition (jTUJl implies = ^ for i — 1, . . . , r. 

Define the set 

B(I) = { f3 = (ft, . . . , fa) | for i = 1, . . . , r, $ is a bijection Si(I) — ► {1, . . . , k} } . 

We have 15(1)1 =h\---l r \ . 

For s = (j, k,l,q) G S'(I) and /5 G S(J), introduce the rational function 

1 a_1 1 

u ^ = _ i n ,(o _ • 

^(s) Z J i=b+l ''ft(s) ''ft-iCs) 

Introduce the rational functions 

Ul = fi V' v E II ^> = E 

lla>6 I'-a.feJ- seS(I) IeP(l,n) 

The function a^ jra defines a rational map 

: C' x C n -> M A [l] , 
called the canonical weight function. 

3.2.2. The canonical weight function was introduced in |SVj . The formula for the 
canonical weight function, presented in Section 13. 2. 1\ is proved in |R,S Vj . 
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3.2.3. Examples. Let n = 2. If Z = (1,1,0,..., 0), then 

. . 1 

£*V(t, z) = — — ^ e 2 ,ie 3<2 VA 1 ® v k2 + — ^ ( (1) e 3A v Al <g> v a 2 

1 

+ TIT) ZTW) 7 e 2,iv Al <8> e 3i2 WA 2 + 7777 7 e 3i2 v Al <8> e 2 ,iWA 2 

(*i ~ z i) 
1 

+ 777 777 7 VAl ® e2 ' ie3 ' 2 ^ + „(2) jiwd : ^ ® • 

(ti - Z2)[t 1 - z 2 ) 

If 1= (2,0,...,0), then 



1 


<*?> 


1 




(tf> 


1 


-z 2 ) 









+ —77 m 1 m m e 2i it) Al (g) e 2 ,iVA 2 

1 2 

+ 777 777) 7 UAl ® e 2,i UA 2 • 

3.3. Hypergeometric solutions to the KZ equations. The master function and 
canonical weight function were introduced in jSVj to solve the KZ equations. The 
hypergeometric solutions to the KZ equations with values in SingM^f/] have the form 



I(z) = [ $(M,A,0 1/K ^,n(M) dt 



where dt = Ajj dt * and 7(2) G C ! x {z} is a horizontal family of /-dimensional cycles 

of the twisted homology defined by the multi- valued function ($>)«, see jS V| IV] . 

In this paper, we will not use the hypergeometric solutions to the KZ equations, 
but the Bethe ansatz for the Gaudin model can be developed studying quasi-classical 
asymptotics of these solutions, |RVj . 



3.4. Bethe ansatz for the Gaudin model. For given I and distinct numbers z±, . . . , z n , 
the problem is to diagonalize simultaneously the Gaudin Hamiltonians, restricted to the 
subspace Sing Ma [i]. 

3.4.1. Theorem ( jRVj . cf. |Bal| lBa2] ). Assume that zeC" has distinct coordinates. 
Assume that t G C' is a critical point of the master function A, I). Then the 

vector uj^ n (t, z) belongs to SingMAfZ] and is an eigenvector of the Gaudin Hamiltonians 
with eigenvalues given by the derivatives of the logarithm of the master function : 

d 

(11) Hi{z) uji tn (t,z) = (— log ${t,z, A, I)) uJi,n{t,z) , i = l,...,n. 
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This theorem was proved in |RVj using the quasi-classical asymptotics of the hyperge- 
ometric solutions of the KZ equations. The theorem also follows directly from Theorem 
6.16.2 in EVl, cf. Theorem 7.2.5 in ISVI. 



3.4.2. Notice that a priori the vector ui tn (t, z) belongs to Ma[/], but if t is a critical 
point, then ui tn (t, z) belongs to SingMA[/]. 

3.4.3. The critical point equations (JUJ) are called the Bethe ansatz equations. 

The values of the canonical weight function at the critical points (with respect to t) 
of the master function are called the Bethe vectors. 

4. Different forms of the Bethe ansatz equations for st 2 and s[ 3 

4.1. Bethe ansatz equations for sfc. The Bethe ansatz equations © for take the 
form 

y. (A s , Ql ) _ y. 2 

W 2^ ,(i)_,(i) ' 

4.1.1. Introduce polynomials 



F(x) = U(x-z s ), G{x) = F(x) £ , ^)=II (*? } 

s=l s=l s i=l 



X 



We have deg F = n, deg G = n — 1, deg P = l x . 

4.1.2. Lemma. Assume that the roots of P are simple. Assume that for any s we have 
z s £ {t x , . . . ,t h }. 

Then the roots of P form a solution to system hltyl if and only if the polynomial 
FP" — GP' is divisible by the polynomial P. 

In other words, the roots of P form a solution to the Bethe ansatz equations Mfy) if 
and only if there exists a polynomial H of degree not greater than n — 2 such that P is 
a solution to the differential equation 

(13) FP" - GP' + HP = . 



The lemma is a classical result due to Stieltjes, see |Stij and Section 6.8 in [Sz] . 
4.1.3. Lemma. There exist unique numbers /ii, . . . , [i n , /ii H — • + /J, n — 0, such that 



H{x) _ \- 1 , \- (Aj,A k ) \ n 
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4.1.4. Corollary. Assume that the roots of P are simple. Assume that for any s 
we have z s {t^ , . . . , } . Then the roots of P form a solution to the Bethe ansatz 
equations fflfy if and only if there exist unique numbers fj,i, . . . , \i n , [i\ + ■ • • + n n = 0, 
such that P is a solution to the differential equation 



(14) P „_ P ,£(^o + P £^ (ft _ E 

1=1 j=l J k^j J 

4.1.5. The corollary provides two ways to describe solutions to the sl 2 Bethe ansatz 
equations. 

The original way: a solution is a collection t = (tf\ . . . , t^) satisfying JEJ. The 
second way: a solution is a polynomial P and numbers /ii, . . . , fi n , n\ + ■ ■ ■ + \i n = 0, 
satisfying JSJ). 

4.1.6. Lemma. For given F,G,H, if at least one of the numbers (Ax, ax), • • • , (A n ,ai) 
is not an integer, then the polynomial solution P to equation \1$ is unique up to mul- 
tiplication by a nonzero number. 

Indeed, if (A s ,ai) is not an integer, then P is the unique (up to multiplication by a 
number) solution, uni- valued in a neighborhood of z s . 

4.2. Bethe ansatz equations for SI3. The Bethe ansatz equations © for 5(3 take the 
form 

h 



\- (A a ,ai) _ 2 1 n • 1 j 

2.^1 ai) .. ^(i) ^(i) r Z-^(i) ^(2) ' 

(15) 



,(1) _ y ^ +0) _ f W ^ _ f (2) 

S=l j = l, j^i fc j fc j jr' = l fc j ^ 



E(A S ,Q!2) 2 1 

,(2) _ 2^ .(2) _ .(2) + .(2) _ .1) ' l = l,...,i2 

S=l j = l, j^i fc j ^jf jr' = l I 'i ""j 

4.2.1. Introduce polynomials F(x) = YH=i ( x ~ z s)> Pi( x ) = Yli=i (4^ — x )> 
= n?=i(*i 2) -^) 5 F i( x ) = P2(x)F(x), F 2 (x) = P 1 {x)F{x\ 

G 1 {x) = Pi(x)F(x) + P 2 (x)F(x)J2^ A ' M|i 



«=i x Zs 



n , , 

G 2 (x) = P;(x)F(x) + />, !./•)/••(,•) ^ 1 



We have deg Fx = n + Z25 deg F 2 = n + /1, deg Gx = n + l 2 — 1, deg C 2 = n + li — 1. 
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4.2.2. Lemma f [MVl] L Assume that the roots tf, tf , t^ t\f of P u P 2 are 

all pair-wise distinct. Assume that for any s, i, we have z s ^ {t\ , . . . , if)}. 

Then the roots of Pi, P 2 form a solution to system < f73)) if and only if for i = 1,2, the 
polynomial F%P" — GiP[ is divisible by the polynomial Pi. 

In other words, the roots of Pi, P 2 form a solution to the Bethe ansatz equations 175)) if 
and only if for i = 1,2, there exists a polynomial Hi of degree not greater than deg Fi — 2 
such that Pi is a solution to the differential equation 

(16) Fi P'i - G t P[ + Hi Pi = . 



Proof. The lemma is a corollary of the same result by Stieltjes, see [StiJ and Section 6.8 
in [Sz]. We sketch the proof. 

Assume that there exist such polynomials Pi, P 2 , Hi, H 2 . Substitute x = tj to equa- 
tion (JH3). Then we get 

Pl\tf) _ Gi{tf) 



mf) mf) ' 

This is exactly the tj-ih equation in Hence the roots of polynomials Pi, P 2 form 
a solution to equations (|15j) . This argument is reversible. □ 

4.2.3. Theorem. Assume that the roots tf, tf, tf\..., t[f of P u P 2 are all pair- 
wise distinct. Assume that for any s,i, we have z s ^ {t[ , ■ ■ ■ ,t\ [ }. Then the roots 
of Pi, P 2 form a solution to system < f73)) if and only if there exist numbers fii, . . . , \i n , 
/ii + ■ ■ ■ + /i n = 0, such that 

(I?) p{'p 2 - p[p> 2 + Pin - p{p 2 jr - PiP' 2 v + 

g—\ S —\ S 

^— ' x — 7 - ^ 7 — 

s=l 

4.2.4. Remark. From (fT7|) one may conclude that 



s=l X Zs k^s Zs Zk 



A*. = y^^lM) + (A,, ai )§(*.) + (A s ,a 2 )^(z s ) = Aio g $(t, z ,A,I)) 

< ' v. — 7, Hi /••'.. ny . 

k^s 

c.f. Theorem 13.4.11 Thus fi s is the eigenvalue of the s-th Gaudin operator at the Bethe 
vector corresponding to the solution t of the Bethe ansatz equations. 

4.2.5. Proof of Theorem \4-2.t\ Let us show that (|T7|) implies (fT5|) . Substitute x = tf 
to (Hlj). Then 

mf) mf) ^(A s ,ai) 



mf) mf) tttf 
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This is the first of equations in Substituting x = tj to (fTTj) we get the second of 
equations in (|T3j) . 

Let us show that (fTljj) imply (fTT|) . Adding the two equations of (fTBj) we get 

FP'(P 2 - FP[P' 2 + FP 1 P^ - FP[P 2 (A ^ Ql) - FP\P'i ( ^ a2) " 

s=l s=l 
- FP[P' 2 + HiPi + P 2 P 2 = . 

Equation ()17|) will be proved if we show that —FP^P^ + H\P\ + P2P2 is divisible by 
PiP 2 - For that it is enough to show that — PP/P^ + H\P\ + P2P2 is divisible by P\ and 
divisible by P 2 . 

From the second of equations in we get 

-FP[P' + H 2 P 2 = -FP X P" + FP 1 Pi V l^fl^i . 

Hence —FP'^P^ + P1P1 + P2P2 is divisible by Pi. Similarly it is divisible by P2. □ 

4.2.6. Lemma [4.2.21 and Theorem 14.2.31 provide three ways to describe solutions to the 
£[3 Bethe ansatz equations. 

The original way: a solution is a collection t = (t^ , . . . , t\ , 4 , . . . , t\ 2 ') satisfying 
(JTHJ) • The second way: a solution is a tuple Pi, P 2 , Hi, H 2 satisfying (fTH|) . The third 
way: a solution is a pair Pi, P 2 , and a set of numbers Hi, . . . , /x n , /ii + • • ■ + /i n = 0, 
satisfying (fTTj) . 

4.2.7. There is a fourth way to describe solutions to the SI3 Bethe ansatz equations. 
Namely, under certain conditions, the Sj-orbits of solutions are in one to one corre- 
spondence with third order linear differential operators with regular singular points 
at Zi, . . . , z n , 00, with prescribed exponents at the singular points, and with a quasi- 
polynomial flag of solutions, see precise statements in [MTV]. Under that correspon- 
dence to a solution t one assigns the differential operator 

D t = (^-ln(^))(^-ln'(^))(A_ ln '(P l)); 
ax P 2 ax Pi ax 

where ln'(/) denotes (df/dx)/f for any /, and T^x) = Yls=i( x ~ z sY Aa ' ai ^ for i = 1,2. 

All singular points of D t are regular and lie in {zi, . . . ,z n , 00}. The exponents of D t 
at z$ are 

0,(A s ,a 1 ) + l,(A s ,a 1 + a 2 ) + 2, 
for s — 1, . . . , n, and the exponents of D t at 00 are 

n n 

—li, —li — (S^ A s — /i«i — / 2 o;2, ai) — 1, —h — (V^ — h®i — hot 2 , ai + a 2 ) — 2. 

s=l s=l 
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The differential equation D t u = has solutions u\,u-2,Uz such that 



ui = Pi, Wr( Ml ,u 2 ) = P 2 Ti , Wr(ui,it 2 ,ii 3 ) = T?T 2 



where Wr (u\, . . . ,Ui) denotes the Wronskian of Ui, . . . , u^. 

We will not use this fourth way to describe solutions to the Bethe ansatz equations 
in this paper. 

5. ON SEPARATION OF VARIABLES FOR sl 2 

In this section we describe Sklyanin's separation of variables for s[ 2 jSkj, following 
exposition in |5tj . 



5.1. Polynomial representation. The space of polynomials Cfx*- 1 -*, . . . , x^\ is identi- 
fied with the tensor product 

M A = M kl ®---®M Kn 

of s[ 2 Verma modules by the linear map 

. . . (x {n) ) jn i-> ei\v Al ® • • • ® e{ n A v An . 

Then the sl 2 action on C[x^\ . . . , x^ n '] is given by the differential operators, 

e « = - x ^d 2 x(i) + (A,, ai )d x(l) , fc« = -2x«^ w + (A,, a,) , /« = *W , 

where c^p) denotes the derivative with respect to a;W. 
The Gaudin Hamiltonians take the form 

Hi(z) = 

-x®x®(d x w -d xU) ) 2 + ((Ai,ai)xW - (A,-, ai )x®)(d xW - d xU) ) + (A,, A,) 



z = 1, . . . , n. 



5.2. Change of variables. Make the change of variables from z%, . . . , z n , x^\ . . . , 

to 2i, . . . , z n , m, t/ 1 ), . . . , y( n_1 ) using the relation 



E 



(,) mzi (t - </*>) 



,_, ' II!' \d - -i ) 

where t is a variable. This relation defines u, y^\ . . . , t/™ -1 ) uniquely up to permutation 
of y^\ . . . , y( n_1 ) unless u = Y^i=i x = 0- The ma P 

(*!, . . . , z n , u, . . . , y^) ^( Zl ,...,z n , x«, . . . , 2 W) 

is an unramified covering on the complement to the union of diagonals = yti) an d 
the hyperplane u = 0. 
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5.2.1. We have 



X W = U 



Us^i ( Z i - Z s) 

dy® 1L ,, (y U) - z s 



dx^ ^ u Ui^(y {j) - y {l) 



and then 



* a !r IL #i (y a) - *«) fl 



.7 = 



5.3. Eigenvectors of the Gaudin Hamiltonians. Assume that we have a common 
eigenfunction of the Gaudin Hamiltonians Hi(z) with eigenvalues /ij. Then the eigen- 
function is annihilated by the operators Hi(z) — /Uj. 

Recall that Hi(z) + • • • + H n (z) = and hence /zi + • • • + /i n = 0. 

Consider the following operators: 

n 1 

^i(*0 = 13 (j) _ ~ - /*<) ' j = 1, • • • ,n - 1 . 

They annihilate a common eigenfunction of the Gaudin Hamiltonians Hi(z) with eigen- 
values Hi. 

5.3.1. Theorem (|Skj. [Si])- variables u,y^\ ... ,y^ n ~ 1 ' , we have 

ts s \ 02 (Aj,ai) 1 (A i? A fc ) 

= + Ej^y— + Ej^i-^ + E — -)• 

1=1 1=1 k^ti 

This is the main point of the separation of variables: the operator Kj(z) depends only 
on and does not depend on other variables y^') and u; this differential operator is 
the same operator for all j; moreover, it is the same operator as in the Bethe ansatz 
equation ([Ti|). 

5.3.2. More precisely, Theorem 15.3.11 claims two identities: 

_ n \ ^ — ^ 

( 19 ) dlu) = Z^/^ TIT) \7 V ~ dx( k )) 2 > 

< 2 °) E s »- = EE _ _ , t) «u - s,<») ■ 

2 — J. 4 — J- fZ-f—% 
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5.4. Canonical weight function. Fix a weight subspace Ma[£] C Ma, I = (h). This 
weight subspace corresponds to the subspace of C[x^\ . . . ,x^] of homogeneous poly- 
nomials of degree l\. The canonical weight function u>i tn becomes the following function 

of T (l) r (") z . , fW 



x (i) 



t {1) -z 

j=l i=l S' Zl 

In variables w, j/W, . . . , y^ n ~ l \ z\, . . . ,z n , , . . . , tj , the canonical weight function is 

< 21 > = Vo..iL) • 

where P(x) = n!Li(4 1 '' ~~ x ) as i n Section EHJ 

5.4.1. Notice that the canonical weight function, as a function of u, y^\ . . . , y^ n ~ l \ is 
the product of functions of one variable. This is another manifestation of separation of 
variables. 

5.4.2. Theorem. Assume that the numbers z = (z±, . . . , z n ) are distinct and the num- 
bers t = (t[ , . . . , t\ ) are distinct. Assume that for any s we have z s {ti, ■ ■ ■ , % }. 
Assume that for such a t, the canonical weight function ty(t, z,u,y), as a function of 
u, y, is an eigenvector of the Gaudin Hamiltonians. Then t is a solution to the Bethe 
ansatz equations Mty) . 

Proof. If 2, u, y) is an eigenfunction, then it is annihilated by the operators Kj(z), 
j — 1, . . . , n — 1. Hence t is a solution of (|12j) by Corollary 14.1.41 □ 

5.4.3. Theorem 13.4. II savs that if t is a solution to the Bethe ansatz equations, then the 
value at t of the canonical weight function is an eigenvector of the Gaudin Hamiltonians 
and is a singular vector. Theorem 15.4.21 gives a converse statement: if the value of the 
canonical weight function at some point t is an eigenvector of the Gaudin Hamiltonians, 
then t is a solution to the Bethe ansatz equations and that value is a singular vector. 

5.5. Lemma. Let at least one of the numbers (A 1; ai), . . . , (A n , «i) be not an integer. 
Let /(z/ (1) ,...,y (n - 1) ) be a polynomial. Assume that for some /ii, . . . , fi n , \i\ + • • • + fi n = 
0, we have Ki(z)f = for i = 1, ... ,n — 1. Then there exists a polynomial P of one 
variable such that 

f(yU,...,y^) = P(yU)...P(y< n -V) . 

Proof. Equation K x {z)f = implies that f{y^\ y^" 1 )) = P(yM) g(y( 2 \ 7/ (n_1) ), 
where P(y^) is a polynomial and g(y ( - suitable function. The polynomial 

P is unique up to multiplication by a number. Applying the same reasoning to g we get 
the lemma. □ 
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6. ON SEPARATION OF VARIABLES FOR S [3 

6.1. Polynomial representation. The space of polynomials C[4^]fc=i' 2 ' 3 of 3n vari- 
ables is identified with the tensor product 

M A = M Al <g> • • • <g> M An 

of 3(3 Verma modules by the linear map 

(4Y (4¥ ■ ■ ■ (4 n) ) J? (4 n V 2 " 



e 2*i e 3^i e 3^2 u Ai ® • • • e^i e 3 " 2 v\ n . 



Then the sl 3 action on C[x^] is given by the differential operators, 

JO - JO - ~(0 1 T Bfl JO - ~(0 

e 2 1 — x 1 , e 3 9 — x 2 -r x 3 u m , e 3 1 — x 



h\ = - 2x\ >d ' (o + 4 J <9 ' (o - 4 ; 9 (o + (A*, «i) , 

x l x 2 x 3 

4 l) = - 2xfd ' (o + _ (0 - 4 l) ^w + ( A *> "2) , 

x 2 x l x 3 

e? 2 = - 4'^w + xfd (i)d (i) - xfd u) - xfd (i) <9 m + (Aj, «i)<9 « , 



4?3 = - 4*^(0 + 4 } 9 « + ( A i; «2)9 _(0 • 

x 2 3 2 

e^ 3 = - 4°^(o - $d wd & + x^d ( i)d 2 {i) - x 2 l) d {i) d (i) - 

' £3 x l ■ l 3 ■•'1 X 2 2 x 3 

- (Aj, a 2 )d u)d w + (A,, «i + a 2 )d (i) . 

Then we have the following formula for the Casimir operator: 

n a,j) = |_ x (0 9 2 +x (*) {d q q ) + ( Ai , ai )0 (0 } x «) + 

+ +*?>(0 «>0«> -9 c,,) -4 j) ^)^) + (A^aOS «,} + 

+ {-4°^w +4°^w + (Aj, a 2 )9 « } {4' ?) + x^d u) } + 
+ {4° +4 ) ^)}{-4 J) ^o) +1?^ + (Aj, a 2 )d x u)} + 
+ {-4°^w - « + ^^w^o - w9 (0 - 

- (Aj, a 2 ) 9,(0 (9,(0 + (Aj, ai + a 2 )9.(o }4 J) + 
+ 4^-4^) - x?d xW d x<J) + xfdjuPp - xfd^d^ - 

>■ a; 3 ^1 ^3 -1-1 a; 2 -1-2 x 3 

- (Aj, a 2 ) <9.co + (Aj, ai + a 2 )d.c») } + 
+{-2x?d x(i) +x^d x(i) -x^d x(i) + (Aj,^)}!-^^) -xfd xU) + (Aj,^t^)}+ 
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+{-2xfd 4) +x ( ? ) d 4) -xfd xf + {Ai,a 3 )}{-x®djp -xfdjp + ( Aj , ftl + 2 ^ 2 ) } . 
Rearranging the terms we get 

(22) = + + + , 

where 

(23) = {K^K^ + ^xfd^-xfxfd^id^-d^ - 
- x?x?{djp - d x o,) 2 - xfx { f(d^ - d^) 2 + 

+ ((A*, «!)xS i} - (Aj, ai )x?) (9 ; ( - 9 .w) + 

+ ((A,, a 2 )4 J) - ( A i> «2)4* ) ) M - ^G)) 



(24) = + (9 - 5 ,,,) (a 4<) - 5 a co) - 

- (/;,. - d a p) 2 + {xfxgdjp - xfxfd xf ) (d xf - djp) + 

+ ((Aj, «i + a 2 )a:3 - (Aj, «i + a 2 )4 ) (9,(0 - 9 w) 



(25) flW) = (xfxf^ - xfxf&p) (9 xli) - d xU) ) + 

+ ((Aj, a 2 )4 J) ^o) - ( A i> a 3 )4 ; 9 x w) - d x u)) , 

(26) n^) = { x f x f - x f x f) (d xf - 9 .o)) . 

The meaning of this decomposition of the Casimir element will be explained in Section 

EH 

6.2. Change of variables. Make the change of variables from z\, . . . , z n , x± , x% , xjj 1 , 

(n) (n) (n) , (1) (1) (1) (n-1) (n-1) (n-1) . 

. . ., x\ ',4 \x\' to 21, . . . ,z„, ui,u 2 ,u 3 , y\ ',y\ J ,yy,...,y\ \y y 2 ',y y 3 using the 
relations 

(i) rrn-i /. _ W)\ 

W E^ = -%^. *= 1.2,3, 

1=1 



* - * nr=i (* - z i) 



where t is a variable. 

Denote a; = {x^\x^\x^\ x^ , xf\ xf ] ), u = (ui,u 2 ,u 3 ), 2/ = (j/i , S/?, 3/3 > 

(n-1) (n-1) (n-l)x 

•••,2/i ,2/2 >2/3 )• 

Relations (|27|) define (u,y) uniquely up to a permutation of y^'s, which preserves 

the lower index, unless Uk = Y^=i = f° r some of fc's. The map 

{z,u,y) -> (2,x) 
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is an unramified covering on the complement to the union of diagonals = y^ and 
the hyperp lanes Uk = 0. 

6.3. Degree on M A [Z]. 

6.3.1. For an index J = ft, ft, ft) and k = 1,2, 3, set J k = j\ H h ft. 

A monomial 

Xj = (x? 5 )^ 1 ^^)^ 

belongs to a weight subspace M A [Z] C M A , I = (Z 1; Z 2 ), if Ji + J3 = Zi and J3 + J2 = ^2- 
We will consider the decomposition 

m a [i) = ©n (/i,/2) M A>d [i) , 

where M A ,d[i] is spanned by all monomials Xj with J\ = l\ — d, J3 = d, J 2 = h — d. We 
say that M A >d [l] consists of elements of degree d. 

6.3.2. In coordinates (u,y), an element belongs to M Ajrf [Z], i.e. has degree d, if it has 
the form 

u i r d uiu i r d f d , 

where f d is a polynomial in y. This polynomial does not depend on y%, . . . , y^ 1 ^ if 
d = 0. An arbitrary element of M A [Z] has degree decomposition: 

TP - ,,h ,.h f /,.(!) ,.(1) „,(«-!) J™" 1 )^ 1 

- u 1 u 2 jo{y l ,y 2 ,...,yi ,y 2 ) + 

min (Zi,/2) 

+ £ «r d «|^-Vd(yf\«5 1 \y?\...,yr i \yr i \yr i) ). 

<2=1 

6.4. Casimir element and degree. The Casimir element is given in ()22 p - ()26| ) . 

Using formulae of Sections 16. 2| consider the Casimir element as an operator acting on 
functions of (it, y) . 

6.4.1. Lemma. 

• The operator Qq'^, given in \2'J\) . preserves the degree introduced in Section 
and does not contain derivatives with respect to variables y^\ . . . , y^ 1 . 

• The operator £1$ , given in \24\j , preserves the degree and annihilates functions 
which do not depend on variables y^, . . . , y^ 1 ^. 

• The operator Vl^q , given in \25\) . increases the degree by one. 

• The operator £2<o , given in \2b\) . decreases the degree by one and annihilates 
functions which do not depend on variables 2/3, . . . , y^ 1 . 

The proof is evident. 
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6.4.2. In the 3n-dimensional space with coordinates (it, y) consider the subspace D 
defined by equations 

(i) (i) 

Vi = Vi , 3 = 1, • • • , n - 1 ■ 
The subspace will be called the main diagonal. 



6.4.3. Lemma. Let F be a polynomial in (u,y) . Apply to F. Then the restriction 
of the function F to the main diagonal equals zero. 

Indeed, the factor (x^x^ — x^x^) in fi<Q is zero on D. 



6.4.4. Let F be an element in M A [l}. Let F = F + F x + . . . and Q^F = (^ (iJ) F) + 
(n' ,|3 ''F)i + . . . be the degree decompositions. 

Lemma. The restrictions to D of the functions (JI^F)q and Qq'^Fq coincide, 

(rt^F) \v = (tij' j) F )\v . 

In other words, the restriction to 2) of the leading term (^ h ^F) can be calculated 
using only the operator Qq applied to F and then restricted to D. 

Proof. The polynomial F does not depend on . . . ,y% . Hence £IqFq = and 
Q^pF = 0. The function Q$pF has degree one. The function (F 2 + F 3 + ...) has 
no degree zero part. The restriction to D of the degree zero part of the function Q^'Fi 
is zero by Lemma f6. 4. 31 This proves the lemma. □ 

6.5. Eigenvectors of the Gaudin Hamiltonians. 

6.5.1. Assume that we have a common eigenfunction F of the Gaudin Hamiltonians 
Hi(z) with eigenvalues /ij, % = 1, . . . ,n. Then the eigenfunction is annihilated by the 
operators H^z) — /ij. Recall that //i + • • • + \i n = 0. 
Consider the following operators: 

- 1 

K M = - m _ - ( H ^ z ) - ' j = i,...,n — i. 

i=i y % 



They annihilate a common eigenfunction of the Gaudin Hamiltonians Hi(z) with eigen- 
values Hi. 
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6.5.2. Decomposition of the Casimir operators into graded components, induces 
the decomposition of the Gaudin Hamiltonians into graded components, 

Hi(z) = Hi(z) + H l {z) + Hi(z) >0 + Hi(z) <0 , 

and the decomposition into graded components 

Kj(z) = Kj(z) + Kj(z) + K 3 (z) >0 + K 3 (z) <0 , 

where 

K j( z )o = Yl (Hi(z)o-m) , Kj(z) = "Tfi ' 

i= i Vi-Zi i= i Vi-Zi 

n j " 1 

k j( z )>o = y^^D — H i( z )>o > ^o(*)<o = -77, — Hi{z) <0 . 

i=i Vx - Zi i=1 y?' - Zi 

6.5.3. Lemma. Lei F = Fo + Ft + . . . be the degree decomposition of the common 
eigenfunction of the Gaudin Hamiltonians, then we have 

(28) (^(z)oF )| s = 0, ./ 1 " 1- 

The lemma is a corollary of Lemma 16.4.41 

These are important equations. Later, under certain conditions, we will show how to 
find Fq from these equations and how to recover F knowing F , see Theorem 16.7.41 

6.5.4. Let F be the common eigenfunction of the Gaudin Hamiltonians, F G Ma[Z], 
I = (h, h)- Then 

where / is a polynomial in y^, y% , j = 1, . . . , n — 1, see Section I6~3l 
Theorem. Equations \2$) have the form 

f7Q , ( gV #7 d*f df A ( A . ttl ) 

1 > V Q ..(i)a..(i) + a..(i) Q „.(i) *..C?) *..(?) o..Cj) 



fyPW flyW V? ~ * 

df (Aj,^) f 1 r (Aj,Afc) 

« (j) 2^ m _ + J 2^ (3) _ v w + 2^ _ 

c/ i/2 i=l i/l z « i=l i/l z « fc^i * K 



. 



Recall that D is defined by the conditions y\ = y 2 for j = 1, . . . , n — 1. 
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6.5.5. As in Sklyanin's Theorem 15.3. 1\ equations ()29j) of Theorem 16.5.41 have three 
interesting properties. 

The first property is that the variables have separated. Namely, the j-th equation 
(Kj(z)oFo)\s) = depends only on variables 2/1,2/2 (which at the end are put being 
equal) and does not depend on other variables , y!f and Mi, u 2 . 

The second property is that the differential operator in ()29|) is the same equation for 
all indices j. 

The third property is that the operator in (|29|) is the same as the operator in the 
Bethe ansatz equation (fTTj) of Theorem 14.2.31 

6.5.6. Proof of Theorem \6.5.4 The theorem follows from Theorem I5.H.11 Indeed, in 
order to prove Theorem 16.5.41 it is enough to prove five identities: 
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i=l k^i 

But the first three identities are corollaries of identity (fTTJJ) . and the last two identities 
are corollaries of identity (|2()j) . □ 

6.6. Canonical weight function. Fix a weight subspace Ma[Z] C Ma, Z = 

In the polynomial representation, the canonical weight function u)i >n becomes a func- 



(i) (i) (i) 
tion m x\ ,x 2 \x 3 , 



(n) (n) (n) 

, Xi , x 2 ,x s , Zx, 



- +(!) f ( 2 ) 



,tj 2 2) . After the 



change of variables of Section IB~2l it becomes a function in ui, u 2 , u$, 2/1,2/2 ,2/3 • • • j 

(n-1) (n-1) (n-1) ,(1) ,(1) .(2) .(2) 

IJl , i/2 ) i/3 J *lj • • • ! Z J15 fc l 5 • • • 1 i l 1 5 fc l 1 • • • l L l 2 • 

We will give a formula for the canonical weight function and its degree decomposition. 
First we prepare notation. 
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6.6.1. For a given d = 1, . . . , min (/i, I2), we will sum certain terms over pairs of or- 
dered subsets (k, m), such that k = (hi, . . . ,kd), rn = (mi, . . . , m d ), where ki, . . . , k d 
are distinct elements of (1, . . . , l\) and m\, . . . , are distinct elements of (1, . . . , I2) 
Such (k, m) will be called d- admissible. The summation over (^-admissible pairs will be 
denoted E d - adni (fc,™) ■ 

To a d-admissible pair (k, m), we assign the following function 

(30) S( fc , TO )(t,«,y) = 

d ^ n— 1 ^(1) n 

= n(,a) _ ,(2) (n r ,(i) _ % )( J) _ jjh) n^s - *«)) • 



6.6.2. Set P l (x) = Y]!Utf ] ~ x), P 2 (x) = nti(f ) " 

6.6.3. Theorem. T/ie canonical weight function and its degree decomposition are given 
by the following formula, 

(31) *(t,z,«,i/) = iffiff x 

ll s =l ^{ZbI^Zb) 



X 



min(Zi,Z2) 

~— ' a! U1U2 ^-^ 

0=1 d— adm(fc,m) 



The theorem is a direct corollary of the definition of the canonical weight function. 

6.7. Comments on Theorem 16.6.31 

6.7.1. Here are the first terms of the degree decomposition 

*(t, 2, 2/) = *(t, z, u, y) + *(t, 2, u, y)i + *(t, 2, n, y) 2 + . . . . 

We have 

w w n;=i Pi(y?)Myi j) ) 

V(t,z,u,y) = u^u'i 



n: =1 pi{z s )p 2 {z s ) 



*(t, z, w, y) x = *(t, z, y) x 

'2 - n-1 (1) (3) 



x i« 2 1? 5 - 1£> (0 rt« _ „co J ft® *•) 



«3_ 

2 fe=l m=l r fc Cm 3=1 i/l A r ™ 2/2 

and so on. 



6.7.2. The degree zero term of the canonical weight function is the product of 
functions of one variable. This is a manifestation of separation of variables. 
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6.7.3. The degree zero term \l/o determines all other terms with d > 0. 

Indeed, knowing \l/ we know the roots . . . , of P\(x) and the roots tf\ . . . , , 
of P2{x). Now for a <i-admissible (k,m), the factor ^o^t^m) has simple combinatorial 
meaning. Namely, we cross out d factors with indices in (k, m) from each Pi, P 2 entering 
\l/o an d then multiply the result by 

d , n-l (1) (j) 

1 tt z k a y3 



n .a) .(2) n 



J-J. £ (1) _J2) 11 ™ / , (1) \ • 

a=l L k a Lm a j=l LLs=l\ L k a Z s) 

Notice that from Pi(y[ ) and i^O/^) we cross out factors depending on y±\ and 
then multiply the result by terms depending on y%\ Hence the larger d, the more 
y% ''-dependent factors participate in \l/ ^(k,m) ■ 

6.7.4. Theorem. Assume that the numbers z = (z%, . . . , z n ) are distinct and the num- 
bers t = (ti , . . . , ti, tf' , . . . ,ti) are distinct. Assume that for any s,i, we have 

z s ^ {ti,...,t\ f^}. Assume that for such a t, the canonical weight function z,u,y), 
as a function ofu, y, is an eigenvector of the Gaudin Hamiltonians. Then t is a solution 
to the Bethe ansatz equations A15\) . 

Proof. Let Hi(z)^(t, z, y, u) = fJ,^(t, z, y, u) for i — 1, . . . , n. 

Since ty(t,z,y,u) is an eigenfunction of the Gaudin Hamiltonians, it is annihilated 
by the corresponding operators Kj(z), j = 1, ...,n — 1. Hence the degree zero term 
(Kj(z)^(t, z, y, u)) is zero. Hence (Kj(z) ^(t, z, y, it))o|x> is zero by Lemma 16.5.31 

Using t, define P\(x), P2(x) as in Section IB. 6. 21 Then the function / = P\{yi ) ^2(2/1 ) 
satisfies equation by Theorem 16.5.41 Then t is a solution to the Bethe ansatz 
equations (jT5j) by Theorem 14.2.31 □ 

6.7.5. Theorem 13.4.11 savs that if t is a solution to the Bethe ansatz equations, then the 
value at t of the canonical weight function is an eigenvector of the Gaudin Hamiltonians 
and is a singular vector. Theorem 16.7.41 gives a converse statement: if the value of the 
canonical weight function at some point t is an eigenvector of the Gaudin Hamiltonians, 
then t is a solution to the Bethe ansatz equations and that value is a singular vector. 

6.8. It would be good to have for g = 5^ an analog of Lemma 1531 



7. Appendix. New form of the Bethe ansatz equations for general 

Kac-Moody algebras 

7.1. Kac-Moody algebras. Let A = (aij)lj =1 be a generalized Cartan matrix, an = 2, 
dij = if and only aji = 0, a^- G Z<o if i ^ j. We assume that A is symmetrizable, 
i.e. there is a diagonal matrix D = diag (d\, . . . , d r ) with positive integers di such that 
B = DA is symmetric. 
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Let g = g (A) be the corresponding complex Kac-Moody Lie algebra (see jKj, §1.2), 
P) C g the Cartan subalgebra. The associated scalar product is non-degenerate on f) * 
and dim 1) = r + 2d, where d is the dimension of the kernel of the Cartan matrix A. 

Let «j G f) *, a( G t) , i = 1, . . . , r, be the sets of simple roots, coroots, respectively. 
We have 

(ai,aj) = didij , (A,a 4 v ) = 2(A, ai)/(a u at) for A G f) * . 

7.2. Bethe ansatz equations for the Gaudin model. Let z = (zi, . . . , z n ) be a 
collection of distinct complex numbers. Let A = (Ai, . . . , A n ), A s G f) *, be a collection 

of -weights and I = (Zi, . . . , l r ) a collection of nonnegative integers. Set I — l\ H \-l r 

and 

t-(t {1) t {1) t (r) t [r) ) 

The Bethe ansatz equations for the Gaudin model associated with this data is the 
following system of algebraic equations with respect to t: 

(32) + rfM y: #4^y = o, 

v ' J-^i ,0) _ /-^ t-^i At) _ As) /L^i Ai) _ M 

s=l L j z s S; k=l L j L k s, s^j L j Ls 

where i = 1, . . . , r, j = 1, . . . , k. 

7.3. The second form of the Bethe ansatz equations. For a function of x we write 
/' = df/dx and ln'(J) = f/f. 

For given z, A, t and i = 1, . . . , r, we introduce polynomials 

3=1 

n 

Fi(x) = {U(x-z s )) J] P 3 (x) , G,(x) 

7.3.1. Lemma ([MVlJ). Assume that the roots t of polynomials P\,...,P r are all 
simple, distinct and different from z\, . . . , z n . 

Then t is a solution of the Bethe ansatz equations \3 b fty if and only if for every i = 
1, . . . , r, the polynomial FiP" — G{P' is divisible by the polynomial Pi. 

In other words, the roots of P\, . . . , P r form a solution of the Bethe ansatz equations 
if and only if for every i = 1, . . . , r, there exists a polynomial Hi of degree not 
greater than deg Fi — 2 such that Pi is a solution to the differential equation 

(33) Fi P" - GiP[ + HiPi = . 



n 

= n^-^) (As * ar> ' 

s=l 

= Fi(x) \n'(T t (x)H P 3 (x)-^^) . 



7 A. The new form of the Bethe ansatz equations. 
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7.4.1. Theorem. Assume that the roots t of polynomials P±,...,P r are all simple, 
distinct and different from z\, . . . , z n . Then t is a solution of the Bethe ansatz equations 
\3ty if and only if there exist numbers /ii, . . . , \x n , \i\ + ■ ■ ■ + \i n = 0, such that 

(34) 5> i5 a,) *± + a,) §5 _ £ (aij a .) 35 + 

" 1 \-(A ; ,A,) . n 



8=1 s k^j 3 K 



This fact may be considered as a generalization of Stieltjes' Lemma 14.1.21 (see also 
lj, and Sec. 6.8 in |SzJ) to an arbitrary Kac-Moody algebra. 



Proof. Let 

n r 

F(x) = H(x-z s ) , P(x) = l[Pi(x) . 

8=1 i=l 

Let us show that (|54*|) implies Multiply (J54*)l by P and substitute for x a root of 
Pi, x = ts. Then 



r il ts J jjVi r 3\ is ) a=l L 8 Z a 

This is exactly the ts -th equation of (JH2J) • 

Let us show the converse. If the zeros of P\, . . . , P r form a solution to (JH2J), then for 
z = 1, . . . , r by Lemma f7. 3. 1| we have 

V ; V 2 P ^p.p 2 IIP, / 



where Hi is some polynomial. 

Add all these equations. Then is proved if we show that 

(36) f^Wi + E(«^«i)p^^ 

i=l i<j 1 3 

is divisible by P. To show that it suffices to show that the expression in ()36|) is divisible 
by Pi, . . . , P r . Let us show that it is divisible by P\. 
Indeed the sum 

H X P X + E ^'-'y^jqv 1 ' 1 ' 

Ki<j % 3 

is divisible by P\. Now for i ^ 1, from the z-th equation in (|33j) we obtain that 

HiP + ( ai , ai )^-PF 

P\P. i 
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is divisible by P\. Hence the expression in is divisible by P±. 

Note that the quotient of (J36j) by P is a polynomial of degree at most n — 2 and 
therefore the quotient of ()36|) by PF has the form 
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